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Abstract. - We introduce and study a novel design for a ratchet potential for soliton excita- 
tions. The potential is implemented by means of an array of point-like (delta) inhomogeneities 
in an otherwise homogeneous potential. We develop a collective coordinate theory that predicts 
that the effective potential acting on the soliton is periodic but asymmetric and gives rise to 
the ratchet effect. Numerical simulations fully confirm this prediction; quantitative agreement 
is reached by an improved version of the theory. Although we specifically show that it is most 
interesting for building Josephson junction ratchets capable to rectify time-symmetric ac forces, 
the proposed mechanism is very general and can appear in many contexts, including biological 
systems. 



Introduction. - Ratchet systems have been the object of intensive studies due to their 
promising applications in biological systems and micro- and nano-technologies 

Recently, a great deal of effort has been devoted to generalize the ratchet mechanism for 
point particles to spatially extended systems 5 . One proposal along this direction relates 
to the existence of net transport in homogeneous extended systems driven by homogeneous 
ac forces [HJEJ. However, although this is a very good and feasible way to induce transport 
with specified properties, it cannot be used to rectify time-symmetric forces, because it relies 
on the breaking of the time symmetry of the ac force. In case the ac force is symmetric, the 
alternative route to rectification is to introduce spatial inhomogeneity. Models in this class 
have been studied theoretically |S] and also implemented in superconducting devices 0^ . In 
this case, a drawback is the difficulty of their fabrication at the micro- or nano-scale, because 
controlling the asymmetry is very complex [3] . An additional factor that has to be considered 
in this context is the interplay between disorder and nonlinearity, which can be of fundamental 
relevance in the design of these new transport devices and particularly when competition of 
scales takes place [TOIITT] . 

In this letter, we present a much simpler design for an extended ratchet that works ir- 
respective of the symmetry of the ac force. Specifically, we focus on the sine-Gordon (sG) 
model, among other reasons because of its important applications to superconducting devices 
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such as long Josephson Junctions (LJJ). In this context, our proposal is based on the inclu- 
sion of point-like inhomogeneities, which correspond to micro-shorts along the LJJ |12U13| . 
Notwithstanding, the mechanism is very general and it can be applied to many other soliton- 
bearing systems where the interaction of kinks with point-like inhomogeneities is similar to 
that occurring in the sG model [14116111 .5] . We begin by introducing the perturbed sG model, 
driven by an ac force in the presence of a periodic array of inhomogeneities. We first study 
analytically this problem by means of a collective coordinate for the motion of the kink center. 
We thus show that under certain conditions the effective kink dynamics is that of a particle 
ratchet. We then study, by numerical simulations of the full system, the behavior of the mean 
velocity as a function of the driving amplitude for different frequencies. Although qualitative 
agreement with the analytical predictions is found, the quantitative comparison is not so good. 
We subsequently improve our theory by introducing the width degree of freedom, whose quan- 
titative success makes clear the physical mechanism needed for a correct description of the 
phenomenon. We conclude by summarizing our main results and pointing out future research 
along these lines. 

Model. - Our starting point is the following perturbed equation: 

<i>tt + 0<t> t - <f> xx + [1 + V(x)] sin0 = Asm(ujt + 5 ), (1) 

where the term V(x) sin(</>) accounts for local inhomogeneities which can be, for instance, 
variations of the critical current in a LJJ, and A sin(ajt + Sq) = f(t) corresponds to an ac bias 
current ^Jj. The kink solution of Eq. 0) in the absence of perturbations, i.e., j3 — V = A = 0, 
is given by the expression 



(x, t) =4 arctan I exp 



x-X(t) 



loVl - v 



(2) 



where X(t) and v = X(t) are the kink position and velocity, respectively, and Iq is a measure 
of the kink width at rest. For the sine-Gordon case Iq = 1, but we choose to leave this as an 
explicit parameter to exhibit the physical relevance of this magnitude; besides, in many other 
soliton-bearing systems collective coordinate equations are the same as in this case but lo ^ 1 
(e.g., the 4> 4 system). The solution in Eq. represents, always in the context of LJJ, a flux 
quantum (fluxon) propagating along the junction. 

The question now is: What is the ideal shape of the function V(x) to turn our system 
into a ratchet device based on the fluxon? It is important to recall that, when driven by 
symmetric ac forces, a damped fluxon can only exhibit oscillatory motion 19 . This problem 
was overcome in previous works, and soliton ratchet behavior has been found for this system 
when the potential of the unperturbed sine-Gordon equation becomes asymmetric |18|. or 
when the system is under an inhomogeneous magnetic field |17j . To our knowledge, none of 
the previous works resorted to inhomogeneities for breaking the symmetry of the system, as 
we now do. Our proposal consists of an array of point-like inhomogeneities, which in the case 
of LJJ can be modelled as delta functions if their length is less than the Josephson penetration 
length. We will choose {x n } to form a periodic, asymmetric array; in particular, for this work 
we have specifically chosen three inhomogeneities per spatial period, 

V(x) = e [S(x — x\ — nL) + 5(x — X2 — nL) + S(x — 23 — nL)] , (3) 



where the parameters should satisfy the constraints a, b, c ~ Iq; a,b < c with a 7^ b, where 
L = a + b + c, a = X2 — Xi, b = X3 — x% and c = L — X3. As will be shown below, for a 
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system satisfying these conditions we obtain net motion with a behavior that resembles very 
closely the one found in ratchet-like systems for point-particles. We have to stress that the 
distances between the delta functions have to be of the same length scale as the kink width; 
otherwise, different behaviors could arise like those demonstrated in |15| . as the interference 
between adjacent deltas is lost. Within that requirement, our proposal is very versatile, as 
it is possible in principle to induce directional motion by using an array whose configuration 
presents more than three inhomogeneities per period, in case a different ratchet potential were 
required. 

Collective Coordinate Theory. - As a first step to justify our choice for the perturbative 
term V(x), we present a simple collective coordinate analysis of its effect on the soliton 
dynamics. The idea of this well-known approximate technique for treating soliton-bearing 
equations is to assume that perturbations affect mostly the motion of the soliton center (and/or 
other parameters, as we will see below) and to proceed to a drastic reduction of the number of 
degrees of freedom by deriving an effective equation for the corresponding collective coordinate 
(see, e.g., ^I] f° r a recent review and further references). One of the easiest procedures to 
derive equations for the collective coordinate is by means of the conservations laws, making use 
of the so-called adiabatic approach, first proposed by McLaughlin and Scott |T2]. Following 
straightforwardly the procedure in this reference, it is a matter of algebra to show, using Eq. 
(J5J as an Ansatz for the perturbed Eq. l(T|l. that the corresponding equation of motion for the 
kink center in the limit of small velocities v 2 <C 1 is 

Qr r 

MoX + pMoX = -—-qAsm((jt + S ), (4) 

where 



2 2 2 1 V ' 

cosh (X — xi — nL) cosh (X — x-i — nL) cosh (X — x% ~ nL)\ 

is the effective potential function, we have for this one time put Iq = 1, and M$ = 8 and 
q = 2tt are the soliton rest mass and topological charge, respectively. 

The potential function given by Eq. © is depicted in Fig. ^i) for the perturbation V{x) 
defined in Eq. ©, with X\ — 0.5, X2 = 1., x% = 2.3 and period L = 4. As can be seen from 
the figure, it corresponds to an asymmetric potential characteristic of ratchet systems; in fact, 
equation is the same as that for a point particle in a rocking ratchet (see |2|). We stress, 
however, that Eq. Q), as it is, describes an inertia! rocking ratchet, similar to those studied 
in |20| . As in this case the dynamics is much more complicated, involving dependencies 
on the initial conditions and on other factors, in the following we restrict ourselves to the 
overdamped case (the common situation in ratchet systems 2 ) by taking j3 = 1. As an 
immediate consequence, within this approximation, we can expect that the soliton center 
should move towards the left. Fig. shows, for one particular value of the frequency, the 
prediction of Eq. (0J in this overdamped approach, confirming this expectation and showing 
the typical window behavior of ratchets (see below). 

Numerical results. - To check the prediction of the simple theory we have summarized 
above, we have carried out numerical simulations of Eq. We have used a Strauss- Vazquez 
numerical scheme |21| with free boundary conditions. The spatial and temporal integration 
steps were Ax = 0.1 and At = 0.01 respectively. The spatial interval for the simulations was 
[—50, 150], with the inhomogeneities arranged periodically according to our three delta unit 
cell in [0, 100]. The system was simulated up to times as long as T — 4000 time units. Finally, 
the numerical representation for delta functions is 1/Ax as usual [THUS*] . 



C/(A) = 2 £ ^ 
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Fig. 1 - a) Effective potential for the kink center within the adiabatic approach, originating in 
the perturbation V(x) defined in Eq. Q with e = 0.8, xi = 0.5, X2 = 1., xs = 2.3 and period 
L = 4. b) Mean kink velocity (dX/dt) vs driving amplitude A for the frequency uj = 0.1. Circles: 
direct numerical simulation of Eq. Q , the line being only a guide to the eye; dashed line: collective 
coordinate approach, Eq. £IJ. 



The simulation results are summarized in Fig. |2 As we may see, the system behaves very 
much like a point particle ratchet (cf., e.g., |22123| ; see also Fig. 5 in 9 ). Indeed, we appreciate 
the existence of windows of motion separated by gaps where the motion is oscillatory with 
mean velocity zero and whose extension increases upon increasing frequency. The explanation 
of these gaps and the observation of "quantized" velocity values as typical signatures of ratchet 
behavior can be found in [21] for the particle-like case and in [H] for the extended system. We 
have also verified that, when the soliton leaves the zone in which the array is contained, its 
motion becomes purely oscillatory, as expected |19) . 

Turning now to a more detailed comparison, we have to admit that the agreement with the 
collective coordinate theory of the previous section is not quite satisfactory. Fig. ^p) makes 
this point clear by showing that neither the number of windows nor, of course, their locations, 
are correctly predicted, even in the simpler low frequency case. Searching for an explanation 
of this problem, we analyzed in depth the simulations, finding out that a possible reason for 
this discrepancy is that the soliton shape changes during its motion along the inhomogeneities 
array (namely, its width is oscillating). This feature can not be accounted for within the 
framework of our theory above and therefore we set out to improve it in the next section. 

Improved Collective Coordinate Theory. - In order to account for the phenomenology 
observed in the simulations, we resort to the generalized travelling wave ansatz for solving our 
problem, first proposed in |25U26| for one and two degrees of freedom. As our starting point 
we rewrite Eq. £Q) in the following way 

= —, 1 p = -—-/34>- s in(4>)V(x) + f(t), (6) 

where ip — (f> and H is the Hamiltonian for the unperturbed problem. Following the procedure 
as in an d assuming that <f> has the form of the so-called Rice Ansatz [27j 

4>{x, t) — 4>k[x — X(t), l(t)] = 4arctan I exp 



x- X(t) 

. W) 
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Fig. 2 - Mean kink velocity (dX/dt) vs driving amplitude A for different frequencies: a) u = 0.015, b) 
uj = 0.05, c) u = 0.1. Other parameters are the same as in Fig.0 Circles: direct numerical simulation 
of Eq. Q, the line being only a guide to the eye; dashed line: improved collective coordinate theory, 
Eqs. 0, 



where l(t) is intended to account for the observed oscillations of the kink width, we find: 

X X XI dU 

M l — + (3M lQ—-M l Q j r = -—-qf{t), (8) 

j j V2 f)TJ 

aMoloj + (3aMoloj + M l — = K mt (l, I, X)-^-, (9) 
where K int (lj,X) = with 

£= ^ + ^«M P + iM (!£ + i). (10, 

M = 8, Iq = 1, a = 7r 2 /12, q — 2-k and U(X, I) has the same form as U(X) in Eq. JSJ) but 
with the denominators recast as cosh 2 (X — Xi ~ nL)/l, i = 1,2, 3. 

In order to make more transparent the physical meaning of these equations, we can change 
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variables by introducing the momentum P(t) — MqIqX /l(t). Our equations become 




(11) 



a[P - 211- 20li] 



I 2 \ P 2 ] 



2l 2 dU 
Mok~dl' 



(12) 




1 + 



We thus see that indeed, in case the kink width oscillates, it necessarily couples to the trans- 
national motion, as the derivative of U(X) in Ea. 111(1 contains I, whereas I is in turn directly 
affected by the momentum. As we observe such oscillations in our simulations, our first 
approach with one degree of freedom must necessarily be incomplete. 

The final step is to compare the improved theory to the simulations. Lacking analytical 
solutions, we have numerically solved the ordinary differential equations JSJl and @, comput- 
ing the kink velocity X(t) and its mean value. As shown by Fig. [2 the comparison between 
our improved collective coordinate theory and the simulations is now excellent, as the win- 
dow numbers are correctly estimated and their locations are very accurately predicted. We 
thus see that although the point particle approximation (collective coordinate X(t)) is enough 
to predict the appearance of a ratchet phenomenon, the correct description of the dynamics 
necessitates one additional degree of freedom, l(t), arising from the fact that the fluxon, the 
"particle" in the ratchet, is an extended object that can show internal oscillations. Strikingly, 
in this way, the interplay of the two degrees of freedom leads eventually to a behavior truly 
indistinguishable from a point-particle ratchet. 

Conclusions. In conclusion, we have proposed and tested in simulations an experi- 
mentally feasible and uncomplicated procedure to build a soliton ratchet using modified long 
Josephson junctions. The main advantage of this system is its simple design, that allows an 
easy implementation by means of indentations of the insulating layer (microshorts). An inter- 
esting feature of our system is its ability to rectify ac forces even if they are time-symmetric, 
something that cannot be accomplished by the homogeneous sG model [HIE]. We have been 
able to show analytically that the physical mechanism responsible for the appearance of the 
ratchet effect is the coupled dynamics of the center and width degrees of freedom, whose 
combined evolution is able to make fluxons (extended objects) behave as point-like particles. 
It is important to remark that the reported phenomenon is robust as we have checked that 
the ratchet effect survives even in the presence of noise [23] . 

Finally, as regards the generality of the procedure presented here, we want to stress that our 
results open the door to many other applications. Indeed, the mechanism for the ratchet effect 
we have found, namely the coupling between translational motion and internal oscillations, will 
be relevant in general for topological solitons, such as those found in the (f> 4 and other nonlinear 
Klein-Gordon models. Such models describe applications in a variety of fields ranging from 
biophysics to pattern forming systems (see, e.g., |11I13( for references). Work along these lines, 
oriented specifically to assess the actual role of this phenomenon related to the macromolecules 
modelled by the <fi 4 equation (see [23 and references therein), is in progress. 
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